Nonequilibrium Work distributions for a trapped Brownian particle in a time 

dependent magnetic field 
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We study the dynamics of a trapped, charged Brownian particle in presence of a time dependent 
magnetic field. We calculate work distributions for different time dependent protocols. In our 
problem thermodynamic work is related to variation of vector potential with time as opposed to 
the earlier studies where the work is related to time variation of the potentials which depends only 
on the coordinates of the particle. Using Jarzynski identity and Crook's equality we show that 
the free energy of the particle is independent of the magnetic field, thus complementing the Bohr- 
van Leeuwen theorem. We also show that our system exhibits a parametric resonance in certain 
parameter space. 

PACS numbers: 05.70.Ln, 05.40.Jc 



Equilibrium statistical mechanics provides us an ele- 
gant framework to explain properties of a broad vari- 
ety of systems in equilibrium. Close to equilibrium the 
linear response formalism is very successful in the form 
of fluctuation-dissipation theorem and Onsager's reci- 
procity relations. But no such universal framework exists 
to study systems driven far away from equilibrium. Need- 
less to say that the most processes that occur in nature 
are far from equilibrium. In recent years there has been 
considerable interest in the nonequilibrium statistical me- 
chanics of small systems. This has led to discovery of sev- 
eral rigorous theorems ,called fluctuation theorems (FT) 
and related equalities 0, S, S, i, i, fl 0, 1, S M, 0]j for 
systems far away from equilibrium. Some of these theo- 
rems have been verified experimentally [HE, [T3L [Tol Hfjj 
on single nanosystems in physical environment where 
fluctuations play a dominent role. We will focus on the 
Jarzynski identity [1] and Crook's equality [j| which deal 
with systems which are initially in thermal equilibrium 
and are driven far away from equilibrium irreversibly. 
Jarzynski identity relates the free energy change(A_F) 
of the system when it is driven out of equilibrium by 
perturbing its Hamiltonian (H\) by an externally con- 
trolled time dependent protocol A(i), to the thermody- 
namic work(W) done on the system, given by 



in 



over a phase space trajectory. Here r is the time through 
which the system is driven. Jarzynski identity is, 
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Crook's equality relates the ratio of the work distribu- 
tions in forward and backward ( time reversed ) paths 
through which the system evolves. This relation is given 



by, 



P f (W) 



= e 0W d 



where, Pf and Py, are the distributions of work along 
forward and backward paths respectively. Here, the dis- 
sipative work — W — W r and W r is the reversible 
work which is same as the free energy difference (A_F) 
between the initial and the final states of the system 
when driven through a reversible, isothermal path. If the 
system is driven reversibly all along the path, the work 
distribution will be 5{W - AF), W d = and P f = P b . 
Thus, the above identities are trivially true for reversibly 
driven system. Jarzynski identity follows from equation 
(3). Crooks relation follows from a more general Crooks 
identity which relates ratio of work probabilities of for- 
ward path and that of the reverse path to the dissipative 
work expended along the forward trajectory. 

In this brief report, we will study the applicability 
of Jarzynski identity and Crooks equality in case of 
velocity dependent as well as time dependent Lorentz 
force which is derivable from a generalised potential, 
U = q(4> — A(t).v). Here, A is time dependent vec- 
tor potential, <f> is scaler potential, q is the charge of 
a particle and v is its velocity. Different time depen- 
dent protocols for magnetic fields are considered. Conse- 
quently, we find that, the free energy difference obtained 
using Jarzynski and Crooks e qual ity complements Bohr- 
van Leeuwen theorem [Tt], HHTHI an d thus we arrive at 
an alternative proof of Bohr-van Leeuwen theorem. This 
theorem states that in case of classical systems the free 
energy is independent of magnetic field and hence the 
thoerem concludes absence of diamagnetism in classical 
thermodyanamical equilibrium systems. We finally also 
show that our system, in presence of ac magnetic field ex- 
hibits parametric resonance in certain parameter regime. 

In an earlier related work [l^, [2(| a charged particle 
dynamics in overdamped limit is studied in the presence 
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of harmonic trap and static magnetic field. The work 
distribution have been obtained analytically for different 
protocols. It is shown that work distribution depends 
explicitly on the magnetic field but not the free energy 
difference (AF). 

The model Hamiltonian for our isolated system is, 



First we have taken magnetic field to vary linearly in 
time, i.e., B = Bq^z. Work distributions for both for- 
ward and backward protocols are obtained. In figure (1) 
we have plotted the distributions Pf(W) and Pj,(—W), 
for forward and backward protocol respectively, which 
are depicted in the insets of figure (1). 
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where k is the stiffness constant of harmonic confinement. 
The magnetic field B(t) is applied in the z direction. The 
x and y components of the vector potential, A x A y are 



given by 



qB(t)y 



and 



qB(t) 



respectively. We have chosen 



symmetric guage here. The particle-environment inter- 
action is modeled via Langevin equation including inertia 
pll |. namely 
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yB(t) + 2yB(t) - kx - Tx + rj x (t), (5) 
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B(t) + 2xB(t) -ky-Ty + ri y {t), (6) 



where T is the friction coefficient and rj x and r\ y are the 
Gaussian white noise along x and y direction respectively. 
This thermal noise has the following properties, 

< rji > = 0; < Vi(t)r]j(t') >= 6 tJ 2Tk B T5(t ~ t'). (7) 

With the above prescription for the thermal noise, the 
system approaches a unique equilibrium state in the ab- 
sence of time dependent potentials. Denoting the pro- 
tocol X(t) = %B(t) the thermodynamic work done by 
external magnetic field on the system upto time r is, 
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(xy — yx)B(i)dt. 



We will like to emphasize that this thermodynamic work 
is related to the time variation of the vector potential 
and can be identified as time variation of magnetic po- 
tential — /i.B, W = — fg n.^-dt, where induced mag- 
netic moment is | (xy — yx) = | (r x v) . To obtain value 
of work and its distribution, we have solved equation (5) 
and (6) numerically using verlet algorithm [22j. We first 
evolve the system upto a large time greater than the typ- 
ical relaxation time so that the system is in equilibrium 
and then apply a time dependent protocol for the mag- 
netic field. We have calculated values of the work for 
10 5 different realisations to get better statistics. The 
values of work obtained for different realisations can be 
viewed as random samples from the probability distri- 
butions P(W). We have fixed friction coefficient, mass, 
charge and ksT to be unity. All the physical parameters 
are taken in dimensionless units. 



FIG. 1: Forward Pf(W) and backward work probability dis- 
tribution Pi(—W) for a ramp magnetic field. 

Using Jarzynski identity (equation (2)) we have com- 
puted free energy difference AF. We have obtained 
(e~^ w ) to be unity (1.0 ± 0.04) implying AF = 0. It 
may be noted that AF = F(B) - F(0), where B is the 
value of the field at the end of the protocol. In the begin- 
ing of the protocol, the value of B is zero. For different 
values of final magnetic field we have obtained AF = 
within our numerical accuracy. This implies that free 
energy itself (and not the free energy difference) is in- 
dependent of the magnetic field, thereby satisfying the 
Bohr- van Leeuwen theorem as stated earlier. We can 
also employ Crook's equality (equation (3)) to determine 
the free energy difference. It follows from Crook's equal- 
ity that Pf and Pb distributions cross at value W = AF. 
This value, where the two distributions cross each other 
(that is, W = 0), can be readily inferred from figure (1). 
This again suggests that, AF = which is consistent 
with the result obtained using Jarzynski identity. 

To strengthen our assertion (that is, the free energy 
being independent of magnetic field) further in figure (2) 
and (3) we have plotted Pf(W) and Pb{— W) for two 
other different protocols as shown in insets of correspond- 
ing figures. For figure (3) we have considered sinusoidally 
varying magnetic field B — B\sinujt in z direction. From 
the crossing point of Pf and Pb we observe that AF = 0, 
consistent with earlier result. 

In figure (4) we have plotted P f (W) and P b {-W)e l3Wd , 
corresponding to the protocol shown in figure (3). Both 
the graphs fall on each other (within numerical error), 
thus verifying Crook's equality. It may be noted that 
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FIG. 2: P f (W)an& P b (-W) for symmetric ramp for B(t). 
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FIG. 4: P/(W) and P b {-W)e l3Wd are plotted together 




FIG. 3: PfiW) and Pb(-W) for oscillatory magnetic field 



reverse protocol also implies reversing the magnetic field 
[23l ] . In all our figures the distribution of work is asym- 
metric and depends on the magnetic field protocol explic- 
itly as opposed to AF. Moreover, all the distributions 
show significant tail in the negative work region. This is 
necessary so as to satisfy Jarzynski identity. 

We now discuss very briefly the occurrence of paramet- 
ric resonance [24| in our system in presence of sinusoidally 
oscillating magnetic field B(t) = B\sinujt. In the param- 



eter range 



qi gl 



Ti > where L\ = 1 



2(fci-rj) 



(see 



4V2LT * " - ^ qt . 

Appendix) , our system exhibits instability. Here k± = — , 
Ti = and q\ = The external parametric mag- 
netic field injects energy into the system and this pump- 
ing is expected to be strongest near twice the systems 
frequency (vLi)- The trajectory of the Brownian parti- 
cle grows exponentially in time also exhibiting oscillatory 
motion at twice the frequency of external magnetic field. 



FIG. 5: x and y Coodinates of the particle and the B(t) are 
plotted as a function of time. 



This is shown in figure (5), where the coordinates of the 
particle and the protocol are plotted as a function of time. 
The parameters are B\ = 60 and u> = 1. For these graphs 
noise strength ksT is taken as one. In presence of this in- 
stability (large variation in coordinate values) it becomes 
difficult to calculate work distributions as it requires large 
number of realisations and better accuracy. Further work 
in analysing the nature of the parametric resonance and 
associated work distributions is in progress. 

In conclusion, by considering the dynamics of a 
trapped charged Brownian particle in a time dependent 
magnetic field we have verified Jarzynski identity and 
Crook's equality. As a by product our result comple- 
ments Bohr-van Leeuwen theorem. Work done on the 
system by external field arises due to the time variation 



4 



of vector potential. This is in contrast to earlier stud- 
ied models where the input energy to the system comes 
from time variation of the coordinate dependent poten- 
tials. Finally, we have discussed very briefly the occur- 
rence of parametric resonance in our system. Our results 
are amenable to experimental verification. 



Redefining T' and q' as Y\ = V /2 and q\ = q' /2 we get, 



£ + [ki - (ri + i qi B lS inujt) 2 }^ = 0. 
Again after transforming t as t 



u = we get, 



9l Sl 



(A5) 



and us as 



I. APPENDIX 

In presence of oscillatory magnetic field B(t) = 
Bisinujt, the mean values of coordinates < x >, < y > 
of the particle (averaged over thermal noise) obey the 
following equation for z =< x > +i < y > 

mz + (r + iqBisinuit)z + (k + i— ^ — cosu>t)z = (Al), 

With k = mki, T = mf, q = mq' the above equation 
becomes 

Q 1 ' B LO 

z + (T' + iq'Bisinujt)z+(ki+i — cosut)z = (A2). 
Now, using the following transformation, 

1 r* 

z(t) = £(i) exp[-- / (r' + iq'B lS inujt)dt}, (A3) 



equation (A2) becomes 

|+ [fci - ~(T' + iq'B lS inujt) 2 }£, = 0. 



(AA) 



o i [L\ + cos2ujiti + iecosuiti]^ = 0, 



(A6) 



where, L\ = 1 - 



2(fei 



, /T ur ■ and e = Irk- For large Sl ' e is 

a small parameter and hence iecostoiti can be treated as 
perturbative term as long as u>\ is far from 2\fL\. The 
condition L\ > 1 should be maintained. Thus ^ can be 
expanded as £ = £o + e £i + ---- Using this in equation(A6), 
we get (keeping only e° order term), 



+ [Li + cos2witi]^ = 0. 



(A7) 



The above equation exhibits parametric resonance 
p^ijwhen u>x m \f~L\. Near resonating frequency, 

Hence, 

0, 



£o goes as ~ e Btl , where s ~ \= . 
z will grow exponentially, if st\ — Y\t > U, i.e. 
fj==(t + ^r) - Tit > O.The condition given above can 



be maintained if ffigjj — Ti > 0. For small amplitude of 
magnetic field, the trajectories of the particles is stable. 
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